Optical Control of the Scattering Length and Effective Range 
for Magnetically Tunable Feshbach Resonances in Ultracold Gases 
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We describe two-field optical techniques to control interactions in Feshbach resonances for two- 
body scattering in ultra-cold gases. These techniques create a molecular dark state in the closed 
channel of a magnetically tunable Feshbach resonance, greatly suppressing optical scattering com- 
pared to single optical field methods. The dark-state method enables control of the effective range, 
by creating narrow features that modify the energy dependence of the scattering phase shift, as 
well as control of the elastic and inelastic parts of the zero-energy s-wave scattering amplitude. 
We determine the scattering length and the effective range from an effective range expansion, by 
calculating the momentum-dependent scattering phase shift from the two-body scattering state. 

PACS numbers: 03.75.Ss 



I. INTRODUCTION 

Ultracold atomic gases with controllable interactions 
arc now widely studied by exploiting collisional (Fesh- 
bach) resonances [IHH , where the interaction strength is 
controlled by means of a bias magnetic field. In a re- 
cent theoretical paper Q , we suggested a general "dark- 
state" optical method for widely controlling the scatter- 
ing length near a magnetic Feshbach resonance, while 
suppressing spontaneous scattering by quantum interfer- 
ence. In this paper, we provide a more detailed treatment 
of the method and show that the molecular dark-state 
method enables control of both the scattering length and 
the effective range in the two-body collisions of ultra-cold 
gases. 

Precision optical control of the scattering length and 
effective range in two-body scattering enables rapid tem- 
poral control and high resolution spatial control of the 
interaction strength near a Feshbach resonance, open- 
ing many new fields of study, such as non-equilibrium 
strongly interacting Fermi gases [f|. For example, the 
natural time scale in a Fermi gas is the "Fermi-time," 
the time tf for an atom at the Fermi surface to move 
a de Broglie wavelength, i.e, Tf = ^f/vf — h/Ep. For 
Fermi energies Ep in the ks x 1 /iK regime, Tf is sev- 
eral /xs. To explore non-equilibrium dynamics on this 
time scale @ requires fast control of interactions, which 
is readily achieved with optical methods. In simulat- 
ing neutron matter, the unitary Fermi gas provides the 
simplest model, where the scattering length is large com- 
pared to the interparticlc spacing. A more realistic model 
of neutron matter can be realized by adjusting both the 
scattering amplitude and the effective range in a trapped 
Fermi gas, to achieve the known ratios for neutrons, 
where the effective range is comparable to the interpar- 
ticle spacing Q- Of great interest, as discussed below, 
is that using optical control, the effective range can even 
be made large and negative. This is especially interesting 
for the narrow Feshbach resonance in 6 Li, where recent 



theory |8| suggests that a Fermi gas can be even more 
strongly interacting than for the broad resonance in the 
unitary regime. 

In contrast to Bose gases, which suffer from three-body 
inelastic processes near a resonance, two-component 
Fermi gas mixtures are stable as a result of the Pauli prin- 
ciple, and can be rapidly cooled to quantum degeneracy 
by evaporation in the resonant regime [l| . Typically, in a 
Feshbach resonance, an external magnetic field controls 
the interaction strength between spin-up and spin-down 
atoms, by tuning the energy of an incoming, colliding 
atom pair into resonance with that of a bound molecu- 
lar state in an energetically closed channel d, [uj ■ Op- 
tical control of Feshbach resonances has been explored 
previously in Bose gases 11 Jl2| and currently is receiv- 
ing substantial attention [13|. Optical Feshbach reso- 
nances (OFR), which employ photoassociation light to 
drive a transition from the continuum of the incoming 
atom pair state to an excited molecular bound state, 
has been proposed and experimentally observed [H, [u3 — 
Il8j . However, light- induced inelastic collisions and the 
accompanying loss limit its practical applicability. Op- 
tical Feshbach resonances for control of higher partial 
waves, such as p-wave scattering of 171 Yb, has been sug- 
gested [l9| and demonstrated recently [2(J ■ OFR also has 
been studied by using a narrow intercombination line of 
a bosonic gas 88 Sr, with the laser frequency tuned far 
away from resonance (2lj . Submicron-scale spatial mod- 
ulation of an inter-atomic interaction has been observed 
in an alkaline-earth atomic condensate (22[. Recently, 
Bauer and coworkers [l3| have used a single optical field 
to control the scattering length near a magnetic Feshbach 
resonance by driving a transition between the resonant 
ground state and an excited molecular state in the closed 
channel. This elegant method enables substantial tuning 
of the scattering length, but a large laser intensity and 
a large frequency detu ning are required for suppressing 
the light-induced loss [13[. The use of electromagnet- 
ically induced transparency (EIT) to control Feshbach 
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resonances was suggested by Harris (23|. Deb [24| has 
suggested that quantum inference between photoassocia- 
tion and magnetic Feshbach s-wave resonance amplitudes 
permits control of the scattering length and suppression 
of inelastic scattering. The method enables control of 
higher order partial waves as well. It is clear that the de- 
velopment of improved quantum interference methods for 
achieving wide tunability of scattering parameters while 
suppressing light-induced loss and heating [f| [24| will 
greatly extend the applicability of optical control meth- 
ods. 




FIG. 1: (color online). Scheme for "dark-state" optical con- 
trol of a Feshbach resonance using two closed-channel molec- 
ular states. Optical fields of frequencies u\ and 0J2 and Rabi 
frequencies Q.\ and Q.2, respectively, couple two singlet ground 
molecular states \g\) and \g2) to the excited singlet molecular 
state |e); Vhf is the hyperfine coupling between the incom- 
ing atomic pair state in the open (triplet) channel and \gi), 
which is responsible for a magnetically controlled Feshbach 
resonance. 



In the paper, we analyze the molecular dark-state 
method proposed in Ref. Q to determine both the op- 
tically controlled s-wave scattering length and the effec- 
tive range from the relative momentum dependence of 
the scattering phase shift. The basic scheme, Fig. [TJ is 
illustrated for a pair of atoms in two hyperfine states 
(denoted spin- up and spin-down) , which undergoes an s- 
wave collision in the ground electronic state triplet molec- 
ular potential (open channel). The hyperfine interaction 
couples the scattering continuum of the open channel to 
a bound singlet vibrational state \g\) in the energetically 
closed channel. An applied bias magnetic field B tunes 
the total energy of the colliding atom pair downward, 
near \g\), producing a collisional (Feshbach) resonance. 
Two optical fields with frequencies oo\ and 002 couple \g\) 
and \g2) to the electronically excited singlet vibrational 
state |e), creating a "dark" state. All three molecular lev- 
els in the closed channel are assumed to have the same 
total nuclear-electron-spin state as \g\), so that both op- 



tical transitions to the level |e) are fully allowed. For 
example, \gz) can be a different singlet vibrational state 
from I <7i ) . To determine the momentum dependence of 
the s-wave scattering phase shift in the presence of the 
light fields, we use a method similar to that employed 
by Fano [2j|, to treat the coupling of an open channel 
continuum to a bound state in an energetically closed 
channel. 

The primary results of this paper show how both 
the s-wave scattering length (eq. I52"j) and the effective 
range (eq. [SZ]) can be controlled using two optical fields, 
which alter the momentum-dependent phase shift of 
the scattering state. This method is applicable to both 
broad and narrow Feshbach resonances, as occur in 
6 Li. In Appendix [Al we present a simple model to 
determine the parameters for Feshbach resonances with 
large background scattering lengths and apply it to 6 Li. 



II. TIME-DEPENDENT SCATTERING STATE 

We will consider first the relevant states for a broad 
Feshbach resonance. As a concrete example, we will use a 
mixture of the two lowest hyperfine states of 6 Li, denoted 
|1) and 1 2). In a bias magnetic field, the atoms interact 
by s-wave scattering, and have a total energy determined 
by the incoming relative kinetic energy, the combined 
hyperfine energies, and the Zeeman energies. 

For a pair of atoms, one each in states |1) and |2), the 
total magnetic quantum number, M = 0, is conserved in 
a bias magnetic field, B z . There are 5 two-atom states 
for M = 0, which can be written in the "interior" singlet- 
triplet basis \S,ms;I,mi). There are two singlet states 
|0, 0; 0, 0), |0, 0; 2, 0), which differ in the total nuclear spin 
I = 0,2, and are degenerate in the absence of hyperfine 
interactions. In addition, there are three triplet states 
1 1, — 1; 1, 1>, |1,0;1,0), and 1 1, 1; 1, 1) . As shown below, 
however, only one singlet state and one triplet state are 
relevant. At high bias magnetic fields, the antisymmetric 
combination of states |1) and \2) is predominantly the 
triplet electronic spin state 1 1 , — 1 ; 1 , 1) . 

A Feshbach resonance in the s-wave channel arises 
when the bias magnetic field tunes the total energy of 
a colliding atom pair in the open triplet channel into res- 
onance with a bound molecular state in the energetically 
closed singlet channel. For 6 Li, the Feshbach resonant 
state, denoted by \gi), is the ground-singlet 38 th molecu- 
lar vibrational state. As this singlet bound state is lower 
in energy than the incoming unbound triplet states at 
zero bias field, the triplet state |1, —1; 1, 1), which tunes 
downward with increasing B-field, is responsible for the 
resonance, which arises from the hyperfine coupling be- 
tween the triplet and singlet channels. The other two 
triplet states tune negligibly or tune upward, and can be 
neglected for the broad Feshbach resonance at 834 G. 

Restricting attention to the 3 interior states, |0, 0; 0, 0), 
|0, 0; 2, 0) , and |1, — 1;1, 1), we consider the origin of the 
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broad and narrow Fcshbach resonances in 6 Li. The hy- 
pcrfine interaction couples each of the singlet states to 
the triplet. However, one combination of |0, 0; 0, 0) and 
|0, 0; 2, 0} has a zero hyperfme matrix element with the 
triplet state [1, — 1; 1, 1) . This uncoupled singlet state 
is |ViV> = (|0, 0; 0, 0) + 2V2|0, 0; 2, 0))/3 and is responsi- 
ble for the narrow Feshbach resonance in 6 Li, near 543 
G [26l [27j . The narrow resonance arises from a second or- 
der coupling of \1pjy) to |1, —1; 1, 1) through the |1, 0; 1, 0) 
triplet state, which is far detuned from the singlet states 
compared to the hyperfme energy. 

To treat the broad Feshbach resonance, we can there- 
fore consider just two states. The first is the combina- 
tion of singlet states which is orthogonal to \iPn), |<?i) = 
(2V2|0,0;0,0) - |0,0;2,0))/3. The second is the triplet 
state |T) = |1, —1; 1, 1), for which the diagonal element of 
the Zeeman-hypcrfine energy is Et = —cihf/2 — B, 
where auF /h = 152.1 MHz is the hyperfme coupling con- 
stant and [Ib is the Bohr magneton, ^s/h ~ 1.4 MHz/G. 

We take the unperturbed two-state Hamiltonian, in the 
absence of optical fields, to be 

Ho = E gi \g 1 ,v 1 )(g 1 ,v 1 \+ E g2 \g 2 ,v 2 )(g2,V2\ 
+ E e \e,v e )(e,v e \ + V HF (\gi)(T\ + \T)( 9l \) 
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where \gi,vi), \g 2 ,v 2 ) arc two ground electronic- 
vibrational states in the singlet molecular potential and 
|e, v e ) is an excited singlet electronic-vibrational state. 
\T, k\) is a triplet continuum state. To model a large 
positive background scattering length a& g , we include in 
the model a near threshold triplet bound state \T, b), 
where E\, = — Ti 2 /{ma 2 g ). We will omit this term for 
small background scattering lengths and for large nega- 
tive background scattering lengths, but it will be impor- 
tant for our discussion of the effective range. In eq.[TJ the 
singlet-triplet hyperfme coupling for the broad 1-2 reso- 
nance in 6 Li is Vhf = (gi\H' HF \T) = -3a HF /(2V3) = 
—hx 131.6 MHz, where H' HF is the hyperfme interaction. 

As the two optical fields used for this method can 
be tuned relatively close to resonance, the isolated reso- 
nance model [lH is adequate, in contrast to single field 
methods employing very large dctunings and high inten- 
sity to avoid optical scattering, where this approximation 
fails [2l|. According to Fig. [T] the two optical fields in- 
troduce a perturbation, 
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H' = -^e-^ t \e,v e )(g 1 ,v 1 \ 



.■^e-^\e,v e ){g 2 M + h.c., (2) 

where we have used the rotating-wave approximation. 
Here, Qi are the Rabi frequencies for the optical transi- 
tions, which include the vibrational Franck- Condon fac- 
tors. 



We take the time-dependent scattering state, for a 
triplet input state of energy E = E F + fi 2 k 2 /m, to be 
of the form, 

= ci \gi, vt) + c 2 \g 2 ,v 2 ) +c e \e,v e ) 
+ c T (k) \T, k) + c T {b)\T, b) + Y, v * k cr(fc') \T, k'), (3) 

where |T, k) is the selected input triplet scattering state 
(in the open channel), with hk the relative momentum 
of the colliding atoms. Note that all of the amplitudes, 
ci, etc., are both fc-dependent and time-dependent. The 
sum denotes the rest of the triplet continuum for k' ^ fc, 
which leads to the principal part of the integral obtained 
in the continuum limit. For r — > 00, the triplet con- 
tinuum background states (box normalization) have the 
asymptotic form 



(r\k) = 



1 sm(kr + 5b g ) 
kr 



V 



(4) 



where Sb g denotes the phase shift arising from the ef- 
fective triplet potential, which corresponds to the back- 
ground scattering length ab g , i.e., 5b g — > — k a^g for k — > 0. 

To obtain the scattering solution, we solve the time- 
dependent Schrodinger equation, 

{H + H')\* E (t))=ih^\* E (t)). (5) 
Taking projections of eq. [5]with (gi,vi\, we obtain 



ci 



E gi a - ig*(k) cr(k) 



1 

h 

■i 9*{ k ) c T{k') - ig* b c T {b) 

k'^k 
O* 

- 1 — — e c e . 
2 



(6) 



where we have defined the hyperfme coupling strengths 



hg(k) = V HF (k\v!) 
%f> = Vhf (b\vi) 



(7) 

(8) 



Here, the overlap integrals are determined from the spa- 
tial wavefunctions (r\v\) and (r\b), which are vibrational 
wavefunctions in the singlet and triplet channel, respec- 
tively, and (r\k) is given by eq. |U Similarly, projecting 
onto (g 2 , v 2 \ yields 



C2 



n 



E„ 2 c 2 



n 



2 AoJzt 



(9) 



while (e, v e \ gives for the excited singlet electronic state 
amplitude 



E e 



7e 



.^1 



ft 2 



ce = -i-^ce-fce + i^-e-^d + i^-e-*"^. (10) 

Here, we have added a phenomenological decay term cor- 
responding to the excited state molecular spontaneous 
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emission rate j e . As the decay rate arises from a spon- 
taneous emission into a broad range of ground singlet 
vibrational states, the arrival rate of molecules into the 
ground states |<7i) and 1 52) is negligible. Finally, project- 
ing onto (T, k"\ and (T, b\ yield 

c T (k' ± k) = -\ (e t + — CT(fc') - ig(k') ci (11) 



and 



c T (b) 



--(E T + E b )c T (b)-ig bCl , (12) 
n 



where we note that the triplet bound state and triplet 
continuum states are orthogonal. 

These time-dependent amplitude equations are easily 
solved. We expect that all of the triplet amplitudes 
will have a time-dependent phase factor cxp(—iEt/h), 
as the asymptotic input and output states arc energy 
cigenstates in the triplet potential. Hence, we take 
cx{k') = br{k') exp(—iEt/h,) for all k' including k and 
c T (b) = b T (b) exp(-iEt/h). Then, eq. HU and E] yield, 

i / b 2 V 2 \ 
b T (k') = -1[Et + — -E) b T (k') 



for k' ^ k and 
b T (b) 



-ig{k')cx e lEt / n 



■(E T + E b -E) b T {b) 



-ig b ci e 



iEt/h 



(13) 



(14) 



It is obvious that we can eliminate the explicit time 
dependent phase factors in eqs. [13] and [14] by taking 
Ci = 61 exp(—iEt/h). Hence, 



b T (k') = --[E T 
and 



2l/2 



h z k 



m 



-E) brW)-ig(U)bi (15) 



b T (b) = -~(E T + E b - E)b T (b) - igbh. (16) 
h 

Using ci = bi cxp{—iEt/K) in eq. [6] wc obtain an equa- 
tion for b\. We eliminate the explicit time-dependent 
phase factor that appears in the b\ equation by taking 
c e = b e cxp[— i(u>i + E/h)t], which yields 



61 



-~(E sl -E)bi-ig*(k)br(k) 
-iJ29*(k')b T (k') 

k'^k 

-ig* b b T (b) + i^b e . 



(17) 



Using c e = b e exp[— i(wi + E/h)t] in eq. [TU] we ob 
tain an equation for b e . We eliminate the explicit time- 
dependent phase factor in that equation by taking c 2 = 
b 2 exp[+i(w 2 — lo\ — E/h)t\. Then, 

n 2 



iA e b e 



le 



(18) 



where, the one-photon detuning A e is given by 

E K — E 

A e = LUl - 



h 



(19) 



Finally, using C2 = b 2 cxp[+i(o;2 — oji — E/Ti)t] in eq. 
we obtain 



b 2 = -18 b 2 + I'-^'be, 

where the two-photon detuning S is given by 

c E — E nr , 

= uj 2 — lo\ - 



h 



(20) 



(21) 



For 6 Li, where the molecular j e ~ 2j spont = 2tt x 11.8 
MHz, and for other atoms with a large spontaneous decay 
rate, we can eliminate the excited state amplitude by 
making an adiabatic approximation, where we assume 
b e << 7 e 6 e /2. In this case, the excited state amplitude 
tracks 61 and b 2 . 



Q,i bi + VL 2 b 2 
2(A e + i 7e /2)" 



(22) 



Wc look for scattering state solutions of the amplitude 
equations where b\,b 2 , br{k) and bx(b) are all zero. Then, 
eq. |2"0"1 yields b 2 = Q,* 2 b e /(28), which with eq . |2"21 yields b e 
in terms of b\ . Then eq. [15] yields 



b T (k' j^k) = 



ftg(fc')6i(fc) 
E - E T -h 2 k ,2 /r 



(23) 



and eq. HHgives 6t(6) = Tig b bi{k)/(E — E T — E b ). Using 
these results in eq. [T7] wc obtain bi(k) in terms of the 
amplitude of the input triplet scattering state amplitude 
b T (k), 



h(k) 



hg*{k) b T (k) 



D(E) 



where 

D(E)=E-E gi -^ E (k) 
Here, the shift is given by 



4[A e + zf + &f] 



(24) 



(25) 



E B (A) = J2 



r\\2 



h'\g(k')\ 



2 I„ 12 



ft \gb 



k'^k 



E -En 



h 2 k>' 



E- E T -E b 



. (26) 



where the sum arises from the coupling of the state \gi) to 
the continuum and the last term arises from the coupling 
of |<7i) to a near threshold triplet bound state. Note that 
we will include this term in the model only when the 
background scattering length a bg is large and positive. 

To convert the sums into integrals, we use Ylk'^k = 
[V/(2tt) 3 ] P J 4:Trk' 2 dk', where P denotes the principal 
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part, and we define the volume-independent coupling 
strength 



hg(k') 



V 



(2nY 



■hg(V) = VHF(%\vi), (27) 



With cq. [30] this yields 



Mr) = b T (k) { (r\k) 



k'^k 



mg(k') g*(k) 
k 2 D(E) 



(r\k>) \ . (32) 



where (k'\vi) is the spatial overlap integral of the vibra- 
tional wavefunction (r\vi) with the continuum normal- 
ized momentum eigenstate. In the limit r — > oo, 



Mk' 



1 



sin(fc'r + SL) 



k'r 



(28) 



To determine the shift, we define E b = —H K. 2 /m, 
where K b = l/a bg for large a bg > with |aj, g | large com- 
pared to the range of the scattering potential. Using 
E = Et + h 2 k 2 /m, and noting that J2 k ,^ k becomes the 
principal part P of the integral over k' , we obtain 



m\g b \ 



p 



(29) 



The asymptotic form of the scattering state, (r\^E(t)), 
in the limit r — > oo, determines the s-wave phase shift. 
At large distance, the singlet and triplet molecular wave- 
functions vanish. Hence, the scattering state is deter- 
mined by the triplet continuum part of the wavefunction, 
(r -> oo|* E (t)) -> ipr(r) \T) exp(-iEt/h), where 

Mr) = b T (k)(r\k) + £ b T (k')(r\k'). (30) 

k'^k 

Using eqs. l23l and 124] we obtain 



From eqs. H [271 and | 

as 



we obtain the scattering state 



Mr ^oc) = MaJ_{sin[fer + 5 bg (k)} 



4-tt m k g* (k) 
Bp) 



oo dk' k' sin[k'r-\-8' b (k')]g(k') 



}. (33) 



To evaluate the principal part appearing in cq. 1331 we 
note from eq. [27] that g(k') is an even function of k' 
since (r\k'), cq. [2H is even in k! . Note that S' bg is an 
odd function of k', i.e., the effective range expansion is 
k' cot[5' bg (k')] = — l/a bg + k' 2 rb g /2 and the right side is 
even in k' . Hence the integrand is even in k' . Further, as 
shown below in the context of our simple model, g{k') has 
only pure imaginary poles arising from the exponential 
form of the molecular bound state of size R. In the con- 
vergent half planes, the sine function yields terms of the 
form exp(ifc'r) — > exp(— r/R) — > 0, which make no con- 
tribution to the integral as r — > oo. Hence, the principal 
part in the limit r — > oo is evaluated as 



1 p roo dk' k' sin[k' r+8' b (k')] g(k') 

2 J-oo k'^-k' 2 

§<?(fc) cos[fcr + 6 bg (k)]. 



(34) 



Using eq.[34]in cq.[33j we obtain finally the asymptotic 
triplet scattering state, 



iPt(t oo) 



b T (k) 1 
sfV kr 



sniffer 



. 2ir 2 mk\g(k)\ 2 _ 
S bg (k)] ^ cos[fcr + 5 bg (k)} 



(35) 



r 



where the input triplet continuum state is of energy E = 
Et + Ti k /m and Et = —a-nFj2 — 2fis B is magnetic 
field dependent. 

To determine the total phase shift A(fc) = A(fc) + 
6 bg (k), we write i()t(t —> oo) = A(k) sin[fcr + A(fc) + 
S bg (k)]/(kr), where A(k) is the resonant part of the 
phase shift and S bg (k) is the background part. Then, 
comparing eq. [35] with A(k) cos A sinffcr + S bg (k)] and 
A(k) sin A cos[fcr + S bg (k)], we obtain 



tanA(fc) = 



2-K 2 mk\g{k){" 
DjE) 



(36) 



where the numerator determines the resonance width, 



which arises from the decay of the dressed molecular state 
into the continuum at a rate T(fc), where 



hT(k) 



2Tr 2 mk\g(k)f 



(37) 



as is readily verified using Fermi's Golden rule. 



III. ZERO ENERGY SCATTERING LENGTH 
AND EFFECTIVE RANGE 

We determine the zero energy scattering length a and 
the effective range r e , from the total phase shift A = 



G 



A + S bg using 



fccot[A(fc)] = -- 



(38) 



With the elementary trigonometric relation between 
cot (A + dbg) and cot A and cot Sb g , we have 



fc cot A 
To expand eq. [38] 



(fc cot A)(fc cot Sbg) - k 2 
k cot Sbg + k cot A 

up to order fc 2 , we use 

D(E) 1 k 2 r e 



fccot[A(fc)] = - 

2-K l vn \g\k)Y 



and 



kcot[S bg (k)} = 

CLbg 



r bg 



(39) 



(40) 



(41) 



with obvious notation. 

After some straightforward algebra, keeping terms up 
to order k 2 (we avoid the zero crossing, assuming 5 + 
dbg 7^ 0), wc obtain the zero-energy scattering length 



a bg 



and the effective range, 



f e a 2 + 2aa bg (a + a bg ) + rb g a 2 
(a + a bg ) 2 



(42) 



(43) 



Note that in the limit a = 0, we have r e = r bg and for 
abg = 0, we have r e = r e as expected. 

We are interested in the zero energy scattering length 
and the effective range near resonance, where the reso- 
nant part of the zero energy scattering length is large 
in magnitude compared to that of the background scat- 
tering length. In the near resonance limit \a\ » \a bg \, 
eq. [43] yields 



2 abg 



2\a bg \ 



1 - 



>a bg 



(44) 



Here, wc ignore the effective range rbg for the background 
scattering states. For small \a,bg\, we can take r e =f e . In 
addition, we show in Appendix [X] using a simple model, 
that when | abg | is large compared to the singlet molecular 
size, the leading contributions to r e from r e exactly cancel 
the ab g -dependent terms in the square bracket. Assuming 
a near resonance condition, where \ab g /a\ is small, we can 
then neglect the correction in the parenthesis. 

In the following, we obtain both the zero-energy scat- 
tering length, the optically-induced inelastic decay rate, 
and the effective range, including the modification arising 
from the optical fields. This is accomplished by expand- 
ing eq.|3n]up to order fc 2 . In this way, we obtain both the 
resonant part of the scattering length a and correspond- 
ing effective range r e . 



For the resonant phase shift, we have using eqs. [40] 
and H5] with E = E T + h 2 k 2 /m, 

E T - E m -Z E (k) + — + D ovt (k) 

kcotA =- 2 JtUi" — ' (45) 

where the optical contribution in the numerator is 



D op t{k) = - 



hl^l 2 



4[A e +*^ + H2f] 



(46) 



Note that both A e and S are dependent on E = Et + 
h 2 k 2 /m and are therefore fc 2 dependent, altering the ef- 
fective range. 

Using fccot A = —1/5 + k 2 r e /2, the resonant part of 
the zero energy scattering length is given by the fc = 
contribution, 



1 = E T -E gi -S £ (0) + D opt (0) 
d 27r 2 m|g(0)| 2 

Since Et = —a^pll — 2fj,s B, we can write 

E T - Eg, - E B (0) = -2fi B {B - B Q 



(47) 



(48) 



where the resonance position includes the shift 

Sb(0). We define the width AB of the resonance in 
terms of the background scattering length, by 



27r 2 m|g(0)| 2 = \a bg \2^ B AB, 



(49) 



where AB is positive by definition. We show in Ap- 
pendix [A] that for a background scattering length that is 
large compared to the molecular size, the energy width 
2fiB AB and the shift 2^(0) are equal in magnitude. 
Eq. [47] determines how the optical fields control the zero 
energy scattering length a = ab g + a. 

In the absence of optical fields, D opt — > 0, the res- 
onant part of the scattering length is then a\B] = 
— \abg\AB/(B — Bqo) and the zero energy scattering 
length takes the usual form [l(| 



i[B] = a bg - \a bg \ 



AB 



B-B n 



(50) 



We see that the zero crossing a[Bo] = occurs at a field 
Bq below (above) resonance for a bg negative (positive). 

The resonant part of the effective range is determined 
from the fc 2 terms in the expansion of eq. 1451 



2 e 



h 2 k 2 3E E (fc) | 2 3r >aptW | 

m K gg) \ k = +K S(gj \ k 
2^ m \g(0)\2 

I _ [-2 AlB (B-B oo )+D opt (0)] 



^ 27r2m|g(0)|" 

Eq. \ST\ can be rewritten as 



f (°) -)- f ' 4- r ° pt 



(51) 



(52) 
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where r e °^ 

-^1,2 



-2(h 2 /m)/(2ir 2 m\g(Q)\ 2 ) arises from the 
h 2 k 2 /m term in eg. 1511 i.e., the relative kinetic energy. 
This term is always present, i.e., even if the shift is inde- 
pendent of k 2 . Using eq. 032 we obtain 



ri 0) 



m[i B AB \a bg Y 



(53) 



which gives the resonant part of the effective range in 
the absence of optical fields when the shift Es(fc) is in- 
dependent of k 2 . We see that broad resonances with 
large background scattering lengths, this contribution to 
the effective range will be small, while it can be large 
for narrow resonances with small background scattering 
lengths [1, HI]. 

The fg term arises from the energy- dependent shift 
Es(fc) and coupling |g(fc)| 2 , which may vary rapidly with 
k 2 when \ab g \ is large, producing large contributions to 
the effective range. Using cq. [47] for the resonant part of 
the scattering length a, this term can be written in the 
form 



k=0 



2n 2 m\g(0)\ 2 a\g{0)\ 2 d(k 2 ) 



(54) 



Finally, the optical fields alter the resonant part of the 
effective range, 



o dD opt (k) | 

Z 9(fc 2 ) l*-'=0 



2tt 2 



™\9( 



(55) 



In Appendix [AJ we determine the optical field inde- 
pendent part of the effective range using a simple model 
for | <Zb g | >> R, where R is the effective size of the sin- 
glet vibrational state To evaluate f' e , eq. HH we use 
cq. IA7l to obtain 



d^ E {k) 



d{k 2 ) 



-2ir 2 m\g(0)\ 2 \ab g \{2e[a bg } ~ 1}. (56) 



fe=0 



and eq. IA4| which gives 
d\g(k)\ 2 



Then, eq. 



d{k 2 ) 
I yields 



KHs(o)|' 



(57) 



k=0 



2\a bg \{2 6[a bg ] - 
2\a bg \ 2 



1} 



2\a bg \' 



— —2a bg — 



(58) 



where the thcta function assures that the first term is just 
—2a bg , for either positive or negative a bg . We see from 
cq. 1581 that r' e exactly cancels the corresponding terms in 
the square bracket of eq. |4"41 

Hence, neglecting the small correction arising from the 
parenthesis in cq. 1441 the effective range for both small 



and large background scattering lengths takes the simple 
form 



f (0) + fopt 
' p I ' p 1 



where r° pt is given by eq. !55l and is given by eq. 1531 In 
the absence of optical fields, wc see that r e 



(59) 



sfl) 



which 



is usually obtained by ignoring the energy dependence of 
the shift and width. 

As D opt (k) in eq.[5S]is a function of E = ET + h 2 k 2 /m, 
we can write 2 dD opt (k)/d(k 2 ) = 2 (fr 2 /m)dD opt (E)/dE. 
Then, since — 2(h 2 /m)/(2Tr 2 m\g(0)\' 



we have 



?(0) 



dD opt {E) 



dE 



k=0 



(60) 



The D opt term enables optical control of the effective 
range, as it can be made to vary rapidly with energy near 
a dark-state resonance. In the following, we systemati- 
cally examine the real and imaginary parts of eq. 1471 and 
eq.[60l 



IV. OPTICAL CONTROL OF THE 
SCATTERING LENGTH 

To find the real and imaginary parts of the zero energy 
(k = 0) scattering length from eq.271 we set a = a bg +a = 
a' + ia" . For this purpose, it is convenient to define the 
magnetic field detuning, 



A = 2 t i B (B-B oo )/n, 

so that -HA Q = -2fi B (B - Boo) = E T - 
Then, we obtain the simple form 

i 

a = a bg - \a bg \f3 



A 



l^il 2 



(61) 

si - MO). 

(62) 



4(5 e + l7c /2) 

where j3 = 2fi B AB/h. Here, we have defined 

IOJ2 

5 P = A e 



AS ' 



(63) 



Note that all the detunings are evaluated for k = 0, i.e., 
with E —> Et = —clhf/2 — 2[i B B the magnetic field 
dependent triplet energy. 

The real part of the scattering length is then given by 

4A rl + |» 1 | 2 r 2( 5 + ( 7e( 5) 2 A 
a - ab ^ KI/? 4(A r 2 + ^ 1 |V4) 2 + ( 7e( 5Ao) 2 ' (64) 

where r 2 = S S e — SA e + | 1 2 /4 has a dimension of fre- 
quency squared. Here, the one-photon detuning, eq. 1191 
is A e = uji — (E e — Et)/?i, while eq. [2~T1 defines the two- 
photon detuning <5 = C1J2 — c^i — (Et — E g2 )/Ti. 
The corresponding imaginary part is 

„__, ,0 7e|»l| 2 £ 2 ,,,, 

a ~ lab9] 2 4(A r 2 + 5|fi 1 | 2 /4) 2 + ( 7e( 5Ao) 2 - (bi)) 
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The imaginary part of the scattering rate causes in- 
elastic loss, which arises from optical scattering, with a 
two-body rate constant i^2(cm 3 /,s) = —8nha"/m in the 
k = limit. 

The dark state method offers many options for 
controlling the scattering length, the inelastic rate, 
the resonance width, and the effective range. These 
include varying the frequencies of the two optical fields, 
choosing the magnetic field detuning and controlling the 
amplitudes of the optical fields in space and time. For 
the initial discussion, we reproduce here the figures from 
our paper 

Fig. [2] shows the real and imaginary parts of the scat- 
tering length as a function of the two-photon detun- 
ing S. We use the parameters for 6 Li: AB = 300 
G, 2(i B /h = 2tt x 2.8 MHz/G, 7e = 2tt x 11.8 MHz, 
and cibg = — 1405 a ; We take f2i = 0.87 e , ^2 = %1e, 
^2 = u eg . 2 , B — B a = 2 G. 




-5 5 



FIG. 2: (color online). Scattering length as a function of the 
effective two-photon detuning 8 in units of 7 e . Real a /dbg 
(Top blue curve) and imaginary a /a bg (Bottom dashed red 
curve). From Ref. Q. 

The dark-state optical control method enables the sup- 
pression of spontaneous scatter, which would cause sub- 
stantial loss and heating if the Oi beam were applied 
alone, as in Ref. Analogous to dark state methods 
for controlling the ratio of absorption to dispersion, we 
can control the ratio a" / a' . Assuming that the resonant 
part of the scattering length a' is large compared to ab g , 
a!' I a! is given by 

a<l 1 7e|»i| 2 <5 2 , . 

a! ~ 2 A [4rl + ( 7e J)2] + ^ |2r 2 J ' {0G) 

where we recall that T 2 = SA e + |ft 2 | 2 / 4 - For A ^ 0, 
we see that loss is suppressed compared to elastic scat- 
tering by the square of the two-photon detuning 5. For 
A = 0, and large Sl 2 , the ratio is — 2^ e (5/ 1 02 1 2 , which can 
be made small for sufficiently large Rabi frequency, Q 2 - 
Fig. El shows the scattering length as a function of il 2 , 
demonstrating the suppression of a" as fl 2 is increased. 
The corresponding ratio a" /a' is shown in Fig. 01 




2 4 



FIG. 3: (color online). Real a ja^g and imaginary a ja^g 
components of the scattering length as a function of 0.2/^^ 
for fli = 5y e , and 8 = 0.05 y e - All other parameters are the 
same as in Fig. [2] The dashed blue line at the bottom is the 
scattering length without the laser fields (magnetic Feshbach 
resonance); The dashed orange line at the top denotes a" — 
0. Inset: Loss ratio between the "dark-state" scheme and a 
typical single laser scheme (where ^2 = 0) as a function of 
£V7 e . From Ref. [|. 




0.2 0.4 



FIG. 4: (color online). The ratio of the imaginary to the real 
part of the scattering length a /a as a function of 0,2 / 7 e for 
fii = 5 7e, and 8 — 0.05 7e . All other parameters are the same 
as in Fig. [3] 



In addition to inelastic loss arising from optical scat- 
tering, there is a small but finite photoassociation rate 
that arises from transitions between the triplet ground 
state and the excited singlet state. The triplet to triplet 
photoassociation transition is far away from resonance. 
For example, choosing the e Li excited singlet vibrational 
state v' = 70, the closest vibrational state for the triplet 
state is v' = 62, which is about 40 GHz away from the 
singlet transition. We will not discuss photoassociation 
here. However, we have shown theoretically that the 
small triplet to singlet photoassociation rate is also sup- 
pressed near the dark-state resonance. 

Finally, the dark-state method produces "artificial" 
narrow Feshbach resonances, which enable rapid changes 
in the scattering length for small changes in the mag- 
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netic field. Fig. [5] shows the results for reasonably large 
value of i7 2 , to clearly separate the broad and narrow 
resonances, which have the usual three-peak structure. 



and 85 /8E = —1/K. Hence, 



J3 



(a) 




FIG. 5: (color online). Scattering length as a function of 
B — Bo for fixed laser parameters tti — 8 7e , ^2 = 12 7e , 
uji = uj egi , ui2 = ^e 9 , • All other parameters are the same as 
in Fig. 2: (a) a ja^g (b) a /a,b g . From Ref. Q. 



V. OPTICAL CONTROL OF THE EFFECTIVE 
RANGE 

In addition to controlling the real and imaginary parts 
of the scattering length, the dark-state method creates 
broad and narrow resonances, as shown in Fig. [5j which 
have different optically controllable effective ranges. The 
dark state method can be applied to both the broad Fes- 
hbach resonance in 6 Li, at 834 G, and the narrow res- 
onance at 543 G. Using both of these, we can explore 
the role of the effective range over a wide range. As dis- 
cussed below in more detail, for fixed two-photon detun- 
ing S — 0, i.e., at the dark-state resonance, the scatter- 
ing length remains at the preselected magnetic field value 
and the inelastic scattering length vanishes, while the ef- 
fective range is dependent on the ratio of the intensities 
of the two optical beams. In this case, the effective range 
can be varied at fixed scattering length, with negligible 
scattering. 

We determine the optically-controlled part of the ef- 
fective range, r e from eq. [6(3 using eq. [46j In this case, 
we note that the one-photon detuning A e defined by 
eq. and the two-photon detuning defined by eq. [5T] 
have a simple ^-dependence, yielding 8A e /8E = 1/H 



d{A e + i le /2 + \n 2 \y(45)}/dE\ k=0 

= (l/h) [l + |fi 2 | 2 /(4<5 2 )] 



which gives 



Nil 



|fi 2 | 2 



4(S e + z 7e /2) 2 



(67) 



where 5 e is defined by eq. [53] with all detunings evaluated 
for k = 0, i.e, E = Et- Taking the real r' e and imaginary 
r" parts of eq. \§7\ we obtain 



r' = f (0) 



|» 1 | 2 (4^ 2 + |» 2 | 2 )[4r|-( 7e( 5) 2 
4[4r 2 + ( 7e< 5) 2 ] 2 



_f(0) 



|n 1 | 2 (45 2 + |fi 2 | 2 )r 27e( 5 



212 



(68) 



(69) 



where T 2 = SS e = 5A e + |ft 2 | 2 /4. 

From eq.l5"3l we recall that re = — h 2 / (m^, B AB\ab g |). 
For the broad resonance in 6 Li at 834 G, where db g = 

-1405 a , AB = 300 G [H, w e have \r^ ] \ ~ la « 
\a,bg\- However, for the narrow Feshbach resonance at 
543 GHH, where AB ~ 0.1 G || and a bg = 62 a , 

ri°^ ~ —7 x 10 4 a , as noted in Ref. [28|. Using the dark- 
state method, the effective range can be widely varied. 

In general, the one- field method |S1 2 | — > will alter the 
effective range as well as the scattering length. Assuming 
2 /4 << A e S, and taking the one-photon detuning 

for simplicity, we find 



|n 2 

A P 



and r' e = f ( e 0) (l 



|fti| 2 / 7 2 ). However, a" £ 0, which can cause substantial 
losses for a one-field method. 

For the two-field method, the imaginary part r' e ' van- 
ishes when either 8 = or T 2 = SA e + |^ 2 | 2 /4 = 0. 
The former corresponds to the narrow peak at the cen- 
ter of Fig. O while the latter corresponds to the two 
side peaks. For the limiting case with il 2 ^ satisfying 
|0 2 | 2 /4 >> SA e , and taking 6 — > 0, we obtain 



..(») 



21 



l« 2 | 2 



(70) 



which shows that the effective range is negative, since 
ri°^ < from eq. 1531 and increases in magnitude when 
the ratio of the Rabi frequencies for the two transitions 
IO1/O2I > 1- Note that when 5 — > for dark-state 
scheme, the loss is negligible and the scattering length 
does not change. The optical fields only modify the ef- 
fective range. 

The effective range as a function of f2 2 is shown in 
Fig. [5] when two-photon detuning is small, <5/ 7e = 0.05. 
It is clear that the effective range can be widely modified 
by the laser fields. 
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FIG. 6: (color online). The real part of the optically con- 
trollable effective range (f' e ) in units of the effective range 
without optical fields (ri 0) ) as a function of Rabi frequency 

0. 2 in units of j e . (i) (Blue oscillatory curve) r'^jrf^ for small 
two-photon detuning 8 = 0.05 j e . (ii) (Red asymptotic curve) 
fe/ri°' for S — > 0. (iii) (Dashed orange curve) r' e /rg = 1, 

1. e., the effective range in the absence of optical fields. 



For the case T 2 = 0, i.e., A e = -|£! 2 | 2 /(4<5), we obtain 
instead, 



..(0) 



1 



l^l^ + l^l 2 ) 

(27e^) 2 



(71) 



In this case, for small two-photon detuning 5 and large 
one photon detuning A e , the effective range is large and 
positive, since ri°^ < 0. 



VI. 



APPLICATION TO 6 LI MAGNETIC 
FESHBACH RESONANCES 



In 6 Li, \gi) is the v — 38 Fcshbach resonance state. We 
can take | (72} to be another lower lying vibrational state, 
which are essentially uncoupled to the triplet state, since 
the nearest state v = 37 is lower in energy than the v = 38 
state by 55.8 GHz, while the v = 36 state is lower by 289 
GHz [30| ■ We employ the bound-to-bound singlet molec- 
ular transition from the ground 1 E+ (N — 0) state to the 
excited A 1 !!^ (N — 1) state. Starting from the v = 38 
Feshbach resonance state, the best Franck-Condon fac- 
tor [3l[ arises in a transition to the v' = 70 vibrational 
state, which we take as |e). The nominal wavelength 
for this transition is 672.66 nm, which is readily accessi- 
ble with a diode laser. A second diode laser excites the 
v = 37 or v = 36 to v' = 70 transition. 

We determine the Rabi frequencies from the known 
dipole transition matrix clement. For the v = 38 — > v' = 
70 transition, the oscillator strength is f eg = 0.035 31 



where f eg = /ip 2 /[3(ea e9 ) 2 ]. Here, n' is the z' compo- 
nent (along the internuclear axis) of the electronic tran- 
sition dipole moment between the selected vibrational 
states, e is the electron charge, and a eg = y/fi/(2m e oj eg ) 
is the electron harmonic oscillator length scale for an 
electronic transition of frequency uj eg . At A = 672.66 
nm, ea eg = 6.9 D (1 D = 1 Dcbyc = 10~ 18 esu-cm.) 
Then fi' = 2.2 D. The corresponding laboratory dipole 

operator is fig — [i' (9, ip), where 9, tp arc the Eulcr 
angles of the molecular internuclear axis with respect to 
the laboratory frame. 

For a AM — 0, ±1 transition in the laboratory 
frame from the N = v = 38 (J = 0) ground 
state to the N = 1 v = 70 (J = 1,M) excited 
state, the transition matrix elements are then all [i = 
Mo/v^ = 1-3 D. The corresponding Rabi frequency is 

ft^Hz) = 4.37 MHz M (D) J /(mW/mm 2 ) , which yields 



Sli = 5.7 MHz y I(mW/mm ), where / is the laser in- 
tensity. This result is nearly identical to that given in 
Ref. [32| for the v' = 68 excited state. 

For the trapped atoms, the applied fields will 
produce an effective light-shift potential, Uls = 
as well as spontaneous scattering at a 



l^l7(4A*a 



ser- 
rate I £l' R 1 2 / (4A 2 aser T S p 0ri t), where rt' R is the Rabi fre- 
quency for the atomic D2 transition. For 6 Li, T spont = 27 
ns and the molecular transition is ~ 1.66 nm (~ —1.1 
THz) red detuned from the free atom transition at 671 
nm. The atomic transition dipole moment is 5.9 D, so 

that fl' R = (5.9/1.3) x fti 
For fix = 10 7 e = 120 MHz, £l R = 540 MHz, the corre- 
sponding free atom scattering rate is only 2.2/s, which is 
negligible for the proposed experiments, where the hold 
time in the weakly interacting regime will << Is. 

For these parameters, the corresponding attractive 
light-shift potential, Uls =66 kHz, or Uls — fcs x 3//K. 
This can be eliminated when necessary using a repulsive 
potential, provided by a spatially-matched blue-detuned 
beam at 532 nm. 



25.8 MHz \J /(mW/mm 



VII. CONCLUSIONS 

The molecular dark-state method permits control of 
two-body scattering parameters, while suppressing light- 
induced inelastic loss and heating compared to single- 
field control techniques. As narrow features are produced 
by these dark-state quantum interference methods, opti- 
cal fields can induce a strong dependence of the scattering 
phase shifts on the relative kinetic energy of the colliding 
atoms, enabling control of the effective range as well as 
the zero-energy scattering amplitude. 
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Appendix A: A Simple Model 

When the background scattering length ab g is large in 
magnitude compared to the molecular size, as it is for 
the broad resonance in 6 Li, it is instructive to evaluate 
the width given by eq. [49] and the shift £ e (k) given by 
eq. 1291 using a simple model. 

We assume that the resonant singlet molecular state 
|t>i) and the triplet bound state \b) have the simple forms 



(r|«i) 
(r\b) 



1 



,-r/R 



V2nR r 



1 



,—r/a b 



y , 2-KCLbg 



(Al) 



where we include the \b) triplet state only for large posi- 
tive scattering lengths, where a\, g » R. Otherwise, this 
state is omitted from the calculation. Hence, there is an 
implied ^[afrg] unit step function in the following. We 
assume that the background scattering states are every- 
where given by eq. 1281 i.e., we ignore the small region of 
rapid oscillation in the deep part of the triplet potential 
well. Then, we easily obtain 



(fcM = 



(1 - a bg /R) 



1 + k 2 R 2 



(A2) 



From eq. IA11 we also have 



(b\vx) = 



2\f dbgR 
R + dbg ' 



(A3) 



where we assume a weakly bound (near threshold) triplet 
state of energy Eb = —h 2 /(maf g ), which arises for large 
positive (background) scattering lengths. Here, we again 
neglect the small region in the deep part of the molecular 
potentials, where the overlap integral of the singlet and 
triplet molecular states oscillates rapidly. For \a,b g \ >> 
R, we then have from eq. [27land lA2[ 



\mk)\ 2 = \v HF \ 2 ^f 



ir 2 1 + (ka bg ) 2 ' 
while eq. [8] and IA3I give for positive a,b g > > R, 



|6»| 2 



abg \abg\ 



(A4) 



(A5) 



Using eq. IA41 the principal part term in cq. [29] is readily 
shown to be — 27r 2 TO|g(fc)| 2 /|a; )£ ,|. Then, we have 



Y m _ m\g b \ 2 \a bg \ 2 „ r „ , 



2ir 2 m\g( 



\a bg \[l + {kabgY 



(A6) 
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TABLE I: Dominant singlet and triplet molecular states in 
the molecular (interior) basis \S M s ; I Mj) for 6 Li Feshbach 
resonances. M is the total magnetic quantum number for a 
pair of colliding atoms. 



Mixture 


M 


Singlet 


Triplet 


1-2 (B) 





(2v/2|00;00) - |00;20))/3 


i i; ii) 


1-2 (N) 





(100; 00) + 2V2|00; 20))/3 


i i; ii) 


1-3 


-1 


100; 2-1) 


i i; io) 


2-3 


-2 


|00; 2 — 2) 


|1 — 1; 1 — 1) 



TABLE II: Triplet energy Et and singlet-triplet coupling 
Vhf for Feshbach resonances in 6 Li. Vhf arises from the 
effective hyperfine interaction o,hf(Ii • Si + I2 • S2), where 
anF/h = 152.1 MHz. Et is the Zeeman- hyperfine energy for 
the given triplet molecular state, with (is the Bohr magneton 
fiB/h — 1.4 MHz/G. For the narrow (N) Feshbach resonance 
in the 1-2 mixture, the coupling is second order in the hy- 
perfine interaction and 1 10; 10} is the dominant off-resonant 
intermediate state. 



Mixture 


Vhf (MHz) 


Et 


1-2 (B) 
1-2 (N) 

1- 3 

2- 3 


-Sa HF /(2V3) = -131.6 
-a HF /(E gi V6) = -5.9 
an f /2 = 76.0 
a HF /V2 = 107.5 


-2/j.b B - a HP /2 
-2u, B B - clhf/2 
— 2/ib B 

-2fl B B + CIHF/2 



where we include a unit step function to indicate that the 
contribution from the triplet bound state is to be used 
only when the background scattering length is large and 
positive. Now, eas. IA4I and IA5I show that m|<7b| 2 |ab 9 | 2 = 
2 x 2ir 2 m\g(0)\ 2 /\a bg \. Hence, 



27r 2 m|g(0)| 2 
\a bg \[l + {ka bg )' 



-{2 9[a bg ]-l}. 



(A7) 



From eq. IA7I and l4"9l we see that for |a& s | >> R, the 
width AB and the shift E_e(0) are related by 



E E (0) = 2(i B AB{2 9[a bg ]-l}. 



(A8) 



Hence, the magnitude of the shift is equal to the twice 
the width (i B AB. 

We can apply this simple model to the three broad 
resonances in 6 Li for 1-2, 1-3, and 2-3 mixtures of the 
three lowest hyperfine states, which are described in Ta- 
ble U and Table [TTJ The measured parameters for the 6 Li 
Feshbach resonances are given in Tabic IIII1 taken from 
ref. [1^] for the broad resonances and from ref. [28| for 
the narrow resonance. Recently, improved 6 Li Feshbach 
resonance parameters have been obtained by using radio- 
frequency spectra of dimer pairs in very low density sam- 
ples, which enables resolution of individual trap-radial- 
vibrational states [33j . 

From eq. ESI for a bg < 0, we have E B (0) = -2(i B AB. 
This result can be used to estimate the energy E gi of 



TABLE III: Feshbach resonance parameters for binary mix- 
tures of the three lowest hyperfine states in 6 Li. The broad 
(B) resonance location -Boo, width AB, and background scat- 
tering length a bg in bohr (do) are taken from ref. [2!|. The 
narrow (N) Feshbach resonance parameters are taken from 
ref. [H. 



Mixture 


Boc(G) 


AB 


abg(a ) 


1-2 (B) 


834 


300 


-1405 


1-2 (N) 


543 


0.1 


+62 


1-3 


690 


122 


-1727 


2-3 


811 


222 


-1490 



TABLE IV: Estimated singlet vibrational energy E gi and size 
R for Feshbach resonances in 6 Li. 



Mixture 


E g Jh (GHz) 


R(a ) 


1-2 (B) 


-2(i B (561G)/h = - 


-1.57 


10.3 


1-3 


-2fi B (568G)//i = 


-1.59 


10.2 


2-3 


-2(i B (562G)/h = - 


-1.57 


10.8 


1-2 (N) 


-2fj, B (570 G)/h = - 


-1.60 





the resonant singlet bound state from the locations and 
widths of the broad resonances given in Table IIIII At 
resonance we have Et(B 00 ) = E gi 



E gi = EriB^) - E B (0) ~ E T {B 



E s (0). Then, 
) + 2(i B AB (A9) 



should be the same for all of the broad resonances. For 
the narrow resonance, we assume that the shift E^ (0) ~ 
0. 

We can also compute the effective size R of the \v\) 
state, eq. IA11 Using eq. |3H1 eq. IA41 and eq. f271 we obtain 



R = 



h 2 (iB AB 
m I a bg 1 1 Vhf I 



(A10) 



Expressing R in bohr units ao, we have, R/clq = 8.4 x 
10 5 AB(G)/[|a 6g (a )||V^(MHz)| 2 ], which also should 
be the same for all of the broad resonances. 

Using the parameters in Table |H] and Table IIIII and 
a HF I '(4/Ub) = 27 G, we obtain from eq. IA9l and eq. IA10I 
the results given in Table IIVI 

The singlet energies obtained from all four resonances 
are nearly identical. Hence, the approximation E^(0) = 
—2(i B AB appears to be reasonably accurate for the large 
negative background scattering lengths in 6 Li, which are 
large in magnitude compared to the size of the resonant 
molecular state. Further, the nearly constant value of 
R validates the scaling given by eq. IA10I based on the 
simple model. 

We can compare the results obtained for R with that 
expected using the overlap integrals for the true wave- 
functions in the singlet and triplet potentials. Using 
ea, IA4l and eq.r^Tl we define an effective size R —> R e ff = 



13 



(tt 2 /\a bg \ 2 )\(k\vi)\l^ . Then, with ip Vl (r) = u v Jr and 
i>Tk(r) = u Tk /r, 



R, 



2n 



eff 



(ka bg ) 2 



dr u vl {r)u Tk (r) 



(All) 



where we take the triplet scattering state to be nor- 
malized so that UTk(f oo) = sin[fc(r — o-bg)] as 
k — > 0. For the simple model, eq. IA11 we have u Vl (r) = 
cxp(— r/R)/( \Z2irR). Taking ut{t) = sin[fc(r — dbg)] 



everywhere and assuming \db g \ » R, we immediately 
obtain R e ff = R as assumed above. We have deter- 
mined the overlap integral of the states obtained for the 
real triplet and singlet potentials (which yield the correct 
highest bound states), using the above normalization for 
ur{r) and the triplet scattering length obtained from the 
triplet scattering state a bg = ar — — 2046 ao. This yields 
yield R e f f ~ 11.4 ao, within 10% of the value R ~ 10.5 ao 
obtained from the broad Fcshbach resonance parameters. 



